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ABSTRACT 

Context. Determining accurate and precise stellar ages is a major problem in astrophysics. These determinations are either obtained 
through empirical relations or model-dependent approaches. Currently, seismic modelling is one of the best ways of building accurate 
stellar models, and therefore providing accurate ages. However, current methods are affected by simplifying assumptions concerning 
stellar mixing processes. In this context, providing new structural indicators which are less model-dependent and more sensitive to 
mixing processes is crucial. 

Aims. We wish to build a new indicator for core conditions (i.e. mixing processes and evolutionary stage) on the main sequence. This 
indicator t„ should be more sensitive to structural differences and applicable to older stars than the indicator t presented in a previous 
paper. We also wish to analyse the importance of the number and type of modes for the inversion, as well as the impact of various 
constraints and levels of accuracy in the forward modelling process that is used to obtain reference models for the inversion. 
Methods. First, we present a method to obtain new structural kernels in the context of asteroseismology. We then use these new 
kernels to build a new indicator of central conditions in stars, denoted t„, and test it for various effects including atomic diffusion, 
various initial helium abundances and various metallicities, following the seismic inversion method presented in our previous paper. 
We then study its accuracy for 7 different pulsation spectra including those of 16CygA and 16CygB and analyse how it depends on 
the reference model by using different constraints and levels of accuracy for its selection 

Results. We observe that the inversion of the new indicator using the SOLA method provides a good diagnostic for additional 
mixing processes in central regions of stars. Its sensitivity allows us to test for diffusive processes and chemical composition mismatch. 
We also observe that using modes of degree 3 can improve the accuracy of the results, as well as using modes of low radial order. 
Moreover, we note that individual frequency combinations should be considered to optimize the accuracy of the results. 

Key words. Stars: interiors - Stars; oscillations - Stars: fundamental parameters - Asteroseismology 


1. Introduction 


Determining stellar ages accurately is crucial when study¬ 
ing stellar evolution, determining properties of exoplanetary 
systems or characterising stellar populations in the galaxy. How¬ 
ever, the absence of a direct observational method to measure 
this quantity makes such determinations rather complicated. 
Age is usually related empirically to the evolutionary stage or 
determined through model-dependent techniques like forward 
asteroseismic modelling of stars. This model-dependence 
is problematic since if a physical process is not taken into 
account during the modelling, we will introduce a bias in the 
age determinations, as well as in the determination of other 
fundamental characteristics like the mass or the radius (see for 
example Eggenberger et al. (2010) for the impact of rotation on 
asteroseismic properties, Miglio et al. (2015) for a discussion 
in the context of ensemble asteroseismology and Brown et al. 
(1994) for a comprehensive study of the relation between 
seismic constraints and stellar model parameters). It is also clear 
that asteroseismology probes the evolutionary stage of stars and 
not directly the age. In other words, we are able to analyse the 
stellar physical conditions but relating these properties to an age 
will, ultimately, always be dependent on the assumptions made 
during the building of the evolutionary sequence of the model. 
A general review of the impact of the hypotheses of stellar 


modelling and of asteroseismic constraints on the determination 
of stellar ages is presented in Lebreton et al. (2014a) and 
Lebreton et al. (2014b). 

In that sense, the question of additional mixing processes 
in the context of seismic modelling is central (Dupret 2008) 
and can only be solved by using less model-dependent seismic 
analysis techniques and new generations of stellar models. 
These new seismic methods should be able to provide relevant 
constraints on the physical conditions in the central regions 
and help with the inclusion of additional mixing in the models. 
In this context, seismic inversion techniques are an interesting 
way to relate structural differences to frequency differences 
and therefore offer a new insight on the physical conditions 
inside observed stars. Erom the observational point of view, 
the high quality of the Kepler and CoRoT data but also the 
selection of the Plato mission (Raueretal. 2014) allows us 
to expect enough observational data to carry out inversions 
of global characteristics. In the context of helioseismology, 
structural inversion techniques have already lead to noteworthy 
successes. They have provided strong constraints on solar 
atomic diffusion (see Basu et al. 1996b), thus confirming the 
work of Elsworth et al. (1990). However, their application in 
asteroseismology is still limited. Inversions for rotation profiles 
have been carried out (see for example Deheuvels et al. 2014, 
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for an application to Kepler subgiants), but as far as structural 
inversions are concerned, one can use either non-linear inversion 
techniques (see Roxburgh 2010, 2015, for an example of the 
differential response technique), or linear inversion techniques 
applied to integrated quantities as is done in Reese et al. (2012) 
and Buldgen et al. (2015). 


ing the adiabatic approximation and spherical symmetry, and ne¬ 
glecting surface integral terms. It is only valid if the stellar mod¬ 
els are sufficiently close to each other. If one is working with the 
structural pair (c^,p), where is the adiabatic squared sound 
velocity and p the density, this relation takes on the following 
form; 


In our previous paper (see Buldgen et al. 2015), we ex¬ 
tended mean density inversions based on the SOLA technique 
(Reese et al. 2012) to inversions of the acoustic radius of the star 
and an indicator of core conditions, denoted t. We also set the 
basis of a general approach to determine custom-made global 
characteristics for an observed star. We showed that applying 
the SOLA inversion technique (Pijpers & Thompson 1994) to 
a carefully selected reference model, obtained via the forward 
modelling technique, could lead to very accurate results. How¬ 
ever, it was then clear that the first age indicator was limited to 
rather young stars and that other indicators should be developed. 
Moreover, the model-dependence of these techniques should be 
carefully studied and there is a need to define a more extended 
theoretical background for these methods. The influence of the 
number but also the type of modes used for a specific inversion 
should be investigated. In the end, one should be able to define 
whether the inversion should be carried out or not, knowing the 
number of observed frequencies and the quality of the reference 
model according to its selection criteria. 

In this study, we wish to offer an answer to these ques¬ 
tions and provide a new indicator for the mixing processes 
and the evolutionary stage of an observed star. We structure 
our study as follows: section 2 introduces a technique to 
obtain equations for new structural kernels in the context of 
asteroseismology and applies it to the («, Ti) and the (m, T) 
kernels, where u is the squared isothermal sound speed, Ti the 
adiabatic gradient and Y the current helium abundance profile. 
Section 3 introduces a new indicator of mixing processes and 
evolutionary stages, which is not restricted to young stars, as 
was the case for the indicator presented in Buldgen et al. (2015). 
Having introduced this new indicator, we test its accuracy using 
different physical effects such as including atomic diffusion 
processes with high velocities (up to 2.0 times the solar micro¬ 
scopic diffusion velocities) in the target, changing the helium 
abundance, changing the metallicity and changing the solar 
mixture of heavy elements. Section 4 analyses the impact of 
the type and number of modes on the inversion results whereas 
Sect. 5 studies how the accuracy depends on the reference 
model. We also tested our method on targets similar to the 
binary system 16CygA and 16CygB using the same modes as 
those observed in Verma et al. (2014) to show that our method 
is indeed applicable to current observational data. Section 6 
summarizes our results and presents some prospects on future 
research for global quantities that could be obtained with the 
SOLA inversion technique. 


2. (u,ri) and (u, F) structural kernels 

2.1. Integral equations for structural couples in the 
asteroseismic context 

It has been demonstrated in Gough & Thompson (1991) that one 
could deduce from the variational principle a linear integral rela¬ 
tion between the perturbations of frequencies and the perturba¬ 
tion of structural variables. This equation is obtained by assum¬ 


V,' Jo p Jo ^ 


dx + — + 0(2), 

Qi 


( 1 ) 


where x — j with R the stellar radius, and where the classi¬ 
cal definition of the relative differences between the target and 
model for any structural quantity s has been used; 


6s(x) _ ^obs(-t:) - ■yref(x) 
^ -^ref 


( 2 ) 


In what follows, we will always use the subscript or superscript 
“obs” when referring to the observed star, “ref” for the reference 
model variables in perturbation definitions and inv for inverted 
results. Other variables, such as the kernel functions, denoted 
without subscripts or superscripts, are of course related to 
the reference model and are known in practice. Finally, one 
should also note that the suffix i is just a index to classify the 
modes. Moreover, since it is clear that some hypotheses are not 
suitable for surface regions, a supplementary function, Q(v) 
was added to model these so-called surface effects. It is defined 
as a linear combination of Legendre polynomials, normalised 
by the factor Qi, which is the mode inertia normalised by the 
inertia of a radial mode interpolated to the same frequency. We 
emphasize here that neither this normalisation coefficient nor 
the treatment of surface effects are uniquely defined and other 
techniques are also used (see for example Dziembowski et al. 
1990; Happen et al. 1991; Basu et al. 1996a). 


The kernels of the couple (p, Ti) were already presented 
in Gough & Thompson (1991) who also mentionned the use of 
another method, defined in Masters (1979) to modify Eq. 1 and 
obtain such relations for the (c^, Ti) couple and also the (A^, c^) 
couple. Other approaches to obtain new structural kernels were 
presented (see for example Elliott 1996; Kosovichev 1999, for 
the application of the adjoint equations method to this problem). 
This latter approach has been used in helioseismology where 
it was assumed that the mass of the observed star is known 
to a sufficient level of accuracy to impose surface boundary 
conditions. In the context of asteroseismology, we cannot make 
this assumption. Nevertheless, the approach defined in Masters 
(1979) allows us to find ordinary differential equations for 
a large number of supplementary structural kernels, without 
assuming a fixed mass *. 

Another question arises in the context of asteroseismol¬ 
ogy: what about the radius? We implicitely define our integral 
equation in non-dimensional variables but how do we relate the 
structural functions, for example c^^,^(r) defined for the observed 
star and c^^f(r) defined for the reference model? What are the 
implications of defining all functions on the same domain in 
X = ^ varying from 0 to 1? It was shown by Basu (2003) 
that an implicit scaling was applied by the inversion in the 
asteroseismic context. The observed target is homologously 
rescaled to the radius of the reference model, while its mean 

' The method of adjoint equations previously described could also be 
used but would require an additional hypothesis to replace the missing 
boundary condition. 
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density is preserved. This means that the oscillation frequencies 
are the same, but other quantities such as the adiabatic sound 
speed c, the squared isothermal sound speed u - j will be 
rescaled. Therefore, when inverted, they are not related to the 
real target but to a scaled one. 

This can be demonstrated with the following simple test. 
We can take two models a few time steps from one another on 
the same evolutionary sequence knowing that they should not be 
that different (here, we consider IM©, main-sequence models). 
We then test the verification of Eq. 1 by plotting the following 
relative difference: 




<$V,- 


V,- 


(3) 


with S' defined as follows: 



The results are plotted in blue in the left panel of Fig. 1, where 
we can see that this equation is not satisfied. However, one could 
think that this inaccuracy is related to the neglected surface terms 
or to non-linear effects. Therefore we carry out the same test 
using the (p, Ti) kernels, plotted in red in the right panel of Fig. 1. 
We see that for these kernels, the equation is satisfied. Moreover, 
when separating the contributions of each structural term, we see 
that the errors arise from the term related to in the first case. 
Using the scaled adiabatic sound speed, however, leads to the 
blue symbols in the right panel of Fig. 1 and we directly see that 
in this case, the integral equation is satisfied. This leads to the 
conclusion that inversion results based on integral equations are 
always related to the scaled target and not the target itself, as 
was concluded by Basu (2003). We will see in Sect. 3 that this 
has strong implications on the structural information given by 
inversion techniques. 


2.2. Differential equation for the (m, Fi) and the (m, Y) 
kernels 

As mentioned in the previous section, the method described in 
Masters (1979) allows us to derive differential equations for 
structural kernels. In what follows, we will apply this method 
to the (M,ri) and the (u,Y) kernels. However, this approach 
can be applied to many other structural pairs such as: (c^, Fi), 
(c^, T), (g,Fi), (g,Y), ..., with the squared adiabatic 

sound speed , g - ^ the local gravity, Fi = '^^e 

adiabatic gradient and Y the local helium abundance. We will 
not describe these kernels since they are straightforward to 
obtain using the same technique as what will be used here 
for the (M,ri) kernels. One should also note that a differential 
equation cannot be obtained for the couple (N^, c^), with the 
Brunt-Vaisala frequency, defined as: ^ . 

without neglecting a supplementary surface term. 

The first step is to assume that if these kernels exist, they 
should satisfy an integral equation of the type given in Eq. 1, 


thereby leading to: 


6vi 

V, 


ri,p r 


Jo M Jo 1 1 


( 5 ) 


From that point, we use the definition of u to express the first 
integral in terms of a density perturbation. This is done using the 
definition of the pressure, F, and the cumulative mass up to a 
radial position, r: 


r* Gmp , 

F = I —::^dr + P^urf 

Jr r 

r 2 

m — I 47rr pdf 

Jo 


(6) 

(7) 


where we will neglect the pressure perturbation at the surface. In 
what follows, we will use the non-dimensional forms P - 
where M is the stellar mass, R the stellar radius and G the gravi- 
tational constant, fh - ^ and p - To avoid any confusion in 
already rather intricated equations, we will drop the hat notation 
in what follows and denote these non-dimensional variables: P, 
m and p. Using Eqs. 6 and 7, one can relate u perturbations to P 
and p perturbations as follows: 


6u 6P 6p 
u P p 


However, using Eq. 6, one can also relate P perturbations to p 
perturbations. Doing this, one should note that the surface pres¬ 
sure perturbation is usually neglected and considered as a so- 
called surface effect. Using non-dimensional variables and com¬ 
bining Eq. 6 and Eq. 7 in Eq. 8, one obtains an expression re¬ 
lating u perturbations solely to p perturbations (of course, this is 
an integral relation due to the definition of the hydrostatic pres¬ 
sure, P). One can use this relation to replace ^ in Eq. 5 and after 
the permutation of the integrals stemming from the definition of 
the hydrostatic pressure perturbation, one obtains the following 
integral equation relating K'^ to K'^ : 



One should be careful when solving this equation since one is 
confronted to multiple integrals, with certain equilibrium vari¬ 
ables associated to x or x. Therefore care should be taken when 
integrating to check the quality of the result. To obtain a differen¬ 
tial equation, we note that it is clear that the equation is satisfied 
if the integrands are equal, meaning that the kernels are related 
as follows: 


. _ m{x)p 


• An}Jc 




r’' K‘ r 

I 

Jo P 


u.r, 

dx 


[f'4 

K‘ r 

1 - ‘ dx 

dx 

Jx X^ 

Jo p 



-K‘. 


w,ri ’ 


( 10 ) 

( 11 ) 
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Fig. 1; Verification of Eq. 1 for a set of 120 modes with the ■ being the 40 radial modes, the x being the 40 dipolar modes and the + 
being the 40 quadrupolar modes. The left plot illustrates this verification for the (p, c^) couple when no scaling is applied to 
whereas the right plot illustrates the same verification for the (p, c^) couple when scaling is applied to c^, as well as the (p, Fi) 

couple where no scaling is needed. 


Given this integral expression, one can simply derive and sim¬ 
plify the expression to obtain a second order ordinary differential 
equation in x as follows: 


(Pk 


-H 


Iny^^p J 

dK 

d^K 

J 

in 

dy (c/y)2 

Iny^^'^p mp 

dK 

in 

' 2yF2p 

dy 

fhp dP 

th dp 

3 dP fhp 

AyP^ dx AyP dx 

4 yF 2 p dx ly^l^P 


K, 


( 12 ) 


where k - k = and y - x . The central boundary 

X~p X~p -c 

condition in terms of k and k is obtained by taking the limit 
of Eq. (12) as y goes to 0. The additional boundary conditions 
are obtained from Eq. 10. Namely, we impose that the solution 
satisfies Eq. 10 at some point of the domain. This system is 
then discretised using a finite difference scheme based on Reese 
(2013), and solved using a direct band-matrix solver. 


Two quality checks can be made to validate our solution, 
the first one being that every kernel satisfies Eq. 10, the second 
one being that they satisfy a frequency-structure relation (As 
Eq. 5) within the same accuracy as the classical structural 
kernels (p, Fi) or (p, c^). We can carry out the same analysis 
as in section 2.1, keeping in mind that the squared isothermal 
sound speed will also be implicitely rescaled by the inversion 
since it is proportional to as is the squared adiabatic sound 
speed, c^. The results of this test are plotted in Fig. 2 as well as 
an example of the verification of the integral equation for the 
kernel associated with the f - 0, n = 15 mode. 


The equation for the («, Y) kernels is identical when using 


the following relation: 



dvi _ 

Vi 

I 

J' K’ypdYdx, 



; du 

Jo 

f K'y^dYdX. 

Jo 

(13) 


meaning that Eq. 12 can simply be transposed using the defini- 

K' „ , K' y 

tions: k — -f- and k - One could also start from Eq. 5, use 

x^p x-p ~ 

the following definition: 


dFi ldln(r{}\ 6P 
Fi \ dln(P) /zxp P 


-H 


/ g/«(ri) \ 

\ dZ jYpp 


5Z, 


pfa(ri) \ 6p ^ p/n(Fi) \ 

I dln{p) Izxp p \ dY ]z,p,p 

(14) 


and neglect the 6Z contribution. This assumption is particularly 
justified if one has spectroscopic constraints on the metallicity. 
Nevertheless, the term associated with 6Z is smaller than the 
three other terms and if one is probing the core regions, the 
^ contribution is already very small. Consequently, all of the 
terms of Eq. 14 are small compared to the integral contribution. 
Still, this assumption is not completely innocent if one wishes 
to probe surface regions. When comparing the ( ^ )yp 


( ain{rp \ 

\ av jz,p,p 


6Y, we notice that their amplitude are comparable 


and that the even often larger. However, we have 

to consider that it will be multiplied by 6Z, which is much 
smaller than dY. Moreover, the functions are somewhat alike in 
central regions and thus there will be an implicit partial damping 
of the dZ term when damping the dY contribution if it is in 
the cross-term of the inversion. One should notice that we can 
control the importance of this assumption by switching from the 
(m, Y) kernels to the (m, Fi) kernels. Indeed, if the error should 
be important, the inversion result would be changed by the 
contribution from the neglected term. In conclusion, in the case 
of the inversion of we will present and use in the next sections, 
this assumption is justified, but this is not certain for inversions 
of helium mass fraction in upper layers, for which only numeri¬ 
cal tests for the chosen indicator will provide a definitive answer. 


Knowing these facts, we can search for the (u, Y) kernels 
and using the previous developments and starting from the 
{p, Y) kernels, directly obtained from (p, Fi) or {p,c^), is more 
straightforward. One should also note that we assume, by using 
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Fig. 2: The left plot illustrates the verification of Eq. 10 for « = 15, ^ = 0 kernel noted where Ju_r, is the right hand side of 
this equation and H is the residual. The right plot illustrates the same test as in Fig. 1 for the (u, Fi) kernels. 
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Fig. 3: Same as Fig. 2 for the (u, Y) kernels. 


Eq. 14, that the equation of state is known for the target. As 
illustrated in Fig. 3, we again test our solutions by plotting the 
errors on the integral equation (Eq. 13), and by seeing how well 
our solution for the £ = 0, n = 15 mode verifies Eq. 10. The 
(m, Fi) and (m, Y) kernels of this particular mode are illustrated 
in figure 4. The kernels associated with u are very similar, 
except for the surface regions where some differences can be 
seen. 

3. Indicator for internal mixing processes and 
evolutionary stage based on the variations of 

u 

3.1. Definition of the target function and iink to the 
evoiutionary effect 

Knowing that it is possible to obtain the helium abundance in 
the integral equation 13, one could be tempted to use it to obtain 
corrections on the helium abundance in the core and therefore 
having insights on the chemical evolution directly. However, 
Fig. 4 reminds us of the hard reality associated with these 
helium kernels. Their intensity is only non-negligible in surface 
regions and it is thus impossible to obtain information on the 
core helium abundance using such kernels. 

Another approach would be to use the squared isothermal 


sound speed, u, to reach our goal. Indeed, we know that 
M = ^ oc and that during the evolution along the main 
sequence, the mean molecular weight will change. Moreover 
the core contraction can also lead to changes in the variation of 
the M profile. Using the same philosophy as for the definition of 
the first age indicator (see Buldgen etal. 2015) and ultimately 
for the use of the small separation as an indicator of the core 
conditions (see Tassoul 1980), we build our indicator using the 
first radial derivative of u. Using u instead of allows us to 
avoid the dependence in Fi which is responsible for the surface 
dependence of To build our indicator, we analyse the effect 
of the evolution on the profile of This effect is illustrated 
in Fig. 5. As can be seen, two lobes tend to develop as the star 
ages. The first problem is that these variations are of opposite 
signs, meaning that if we integrate through both lobes, the 
sensitivity will be greatly reduced. Therefore, we choose to base 
our indicator on the squared first radial derivative: 



where f(r) is an appropriate weight function. First, we will con¬ 
sider that the observed star and the reference model have the 
same radius. The target function for this indicator can easily be 
obtained. We perturb the equation for f„ and use an integration 
by parts to relate the perturbations of the indicator to structural 
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Position x=r/R Position x=r/R 


Fig. 4: Structural kernels for the n — 15, f = 0 mode associated with the (m, Fi) couple on the left hand side and with the (m, Y) 

couple on the right hand side. 


perturbations of u: 



The last term on the second line is not suitable for SOLA inver¬ 
sions, given the neglect of surface terms in the kernels. Hence we 
define the function / so that /(O) = f(R) = 0, thereby cancelling 
this term. This leads to the following expression; 



where 




-2m d 
iu dr 



( 18 ) 


The weight function /(r) must be chosen according to a number 
of criteria: it has to be sensitive to the core regions where the 
profile changes; it has to be of low amplitude at the boundaries 
of the domain, allowing us to do the integration by parts neces¬ 
sary to obtain 6u in the expression; it should be possible to fit the 
target function associated with this f(r) using structural kernels 
from a restricted number of frequencies; moreover, Eq. (18) 
being related to linear perturbations, it is clear that non-linear 
effects should not dominate the changes of this indicator^. We 
also know that the amplitude of the structural kernels is 0 in the 
centre, so /(r) should also satisfy this condition. 

^ Otherwise, using the SOLA technique, which is linear, would be im¬ 
possible. 


We define the weight function as follows: 

/(r) = r“(r-/?)-'exp|-r(^^^^)"j, (19) 

which means that we have 4 parameters to adjust. The case of 
ropt is quickly treated. Since we know that the changes will be 
localised in the lobes developing in the core regions, we chose to 
put ropt = 0. a and A should be at least 1, so that the integration 
by part is exact and the central limit for the target function and 
the structural kernels is the same, y depends on the effects of 
the non-linearities. However, since we will have to perform a 
second derivative of m a more practical concern appears: we 
don’t want to be influenced by the effects of the discontinuity 
at the boundary of the convective envelope. At the end of the 
day, we use the following set of parameters: a - 1, A - 2, 
y -1 and rgpt - 0. One could argue that the optimal choice for 
ropt would be at the maximum of the second lobe or between 
both lobes to obtain the maximal sensitivity in the structural 
variations. These values were also tested, but the results were a 
little less accurate than using ropt - 0 and they involved higher 
inversion coefficients and therefore higher error magnification. 
We illustrate the weighted profile obtained for this optimal set of 
parameters in the right-hand side panel of Fig. 5. Furthermore, 
Eq. 17 is satisfied up to 5%, so we can try to carry out inversions 
for this indicator. It is also important to note that for the sake 
of simplicity, we do not chose to change the values of these 
parameters with the model, which would only bring additional 
complexity to the problem. 

Thanks to the target function Tt^, we can now carry out 
inversions for the integrated quantity using the linear SOLA 
inversion technique (Pijpers & Thompson 1994). Eirst, let us 
take some time to recall the purpose of inversions and our 
adaptation of the SOLA technique to integrated quantities. 
Historically, inversions have been used to obtain seismically 
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constrained structural profiles (Basu et al. 1996b) as well as ro¬ 
tation profiles (Schou et al. 1994) in helioseistnology. However, 
all these methods are not well suited for the inversions we wish 
to carry here. As discussed in Buldgen et al. (2015), the SOLA 
inversion technique, which uses a “kernel matching” approach 
is well suited to our purpose. Indeed, this approach allows us to 
define custom-made target functions that will be used to build a 
cost function, here denoted ST. In the case of the f„ quantity, one 
has the following definition; 




[/fAvg-TlJ dx+jBj + tan(0) ^(ao-i) 

N 

J]ci-k 


( 20 ) 


where is the so-called averaging kernel and f^cross is the so- 
called cross-term kernel defined as follows for the (m, Li) struc¬ 
tural pair (For (u, Y), one would replace K‘^ ^ with K‘^ y ^nd Ky ^ 

with 4 „): 


N 

f^Avg - 

i 

N 

= ( 22 ) 

i 

9 and p are free parameters of the inversion. 6 is related to the 
compromise between the amplification of the observational error 
bars (cTi) and the fit of the kernels, whereas p is allowed to vary 
to give more weight to the elimination of the cross-term kernel. 
In this expression, N is the number of observed frequencies, the 
Ci are the inversion coefficients, which will be used to determine 
the correction to be applied on the value. 77 is a Lagrange 
multiplier and the last term appearing in the expression of 
the cost-function is a supplementary constraint applied to the 
inversion which is presented in Sect. 3.2. 


If the observed target and the reference model have the 
same radius, the inversion will measure the value of f„ for 
the observed target. However, if this condition is not met, the 
inversion will produce a scaled value of this indicator. By 
defining integral equations such as Eq. 15, or even Eq. 1, we 
have seen in Sect. 2.1 that we made the hypothesis that both 
target and reference model had the same radius. However, as the 

frequencies scale with , the inversion will preserve the 

mean density of the observed target. Therefore, we are implicitly 
carrying out the inversion for a scaled target homologous to the 
observed target, which has the radius of the reference model but 
the mean density of the observed target. Simple reasoning shows 

US that the mass of this scaled target is: Mtar - Thus, 

as scales as M^, there is a difference between the target value 
and the measured value, f"". We can write the following 
equations: 


24 as a criterion to determine whether the inversion was success¬ 
ful or not. 


3.2. Non-linear generalisation 


We present in this section a general approach to the non¬ 
linear generalisation presented in Reese et al. (2012) and 
Buldgen etal. (2015) for any type of global characteristic that 
can scale with the mass of the star. We can say that the frequen¬ 
cies scale as Of course, they scale as the mean density, 

namely However, since the inversion works with a fixed 
radius and implicitely scales the target to the same radius as the 
reference model, the iterative process associated with this non¬ 
linear generalisation will never change the model radius. There¬ 
fore, we do not take this dependence into account and simply 
work on the mass dependence. Let A be a global characteristic, 
related to the mass of the star. It is always possible to define a 
factor k so that: 

A oc /. (25) 


And we obtain directly: 
— -k— 


(26) 


However, using the definition of the inverted correction of A, one 
has: 



(27) 


with the Ci the inversion coefficients. Using the same reasoning 
as in Buldgen et al. (2015) and Reese et al. (2012), we define the 
inversion as “unbiased” (term which should not be taken in the 
statistical sense.) if it satisfies the condition; 

(28) 


Now we can define an iterative process, using the scale factor 
Sa = we will scale the reference model (in other words, 

multiply its mass and density by s^, its pressure by s'*, leaving E 1 
unchanged) and carry out a second inversion, then define a new 
scale factor and so on until no further correction is made by the 
inversion process. In other words, we search for the fixed point 
of the equation of the scale factor. At the /* iteration, we obtain 
the following equation for the inversion value of A: 


-4inv ~ -^ref-^ 


.k-l 


+ + k(l - Sj) 


(29) 


This can be written in terms of the scale factor only, noting that 
at the iteration, 4 ^ = s^,,: 

■' ’ 7+1 


f' 


fObs 


m 2 

K. 

ref 


tar 

.inv 


■2 ’ 
tar 
.obs 


M: 

f 


R,. 


(23) 

(24) 


where we have used the definition of Mtar to express the mass 
dependencies as radius dependencies. Therefore, we will use Eq. 


^j+i 


„k-l 


Y,Ci^ + - Sj) 


The fixed point is then easily obtained: 

■^opt “ "f" 1 ; 


(30) 


(31) 
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Fig. 5: Left panel: structural changes in the scaled ^ profile with the evolution. The models are main sequence I.OMq models. 
Right panel: the same effects as seen when is weighted according to Eq. 19. 


and can be used directly to obtain the optimal value of the in¬ 
dicator A. We can also carry out a general analysis of the error 
bars treating the observed frequencies and Ainv as stochastic vari¬ 
ables: 


t'obs,/ — t^ob.s./(l "i" (32) 

^inv = ^inv(l + ^a), (33) 


with e, and being the stochastic contributions to the variables. 
Using the hypothesis that e, « 1 we easily obtain the following 
equation for the error bars: 


— ^ref 


1 

k 


yObs 




k-\ 



(34) 


We note that in the particular case of the indicator k - A. 
Indeed, t,, oc Afi whereas v oc 


3.3. Tests using various physicai effects 

To test the accuracy of the SOLA technique applied to the 
indicator, we carried out the same test as in Buldgen et al. 
(2015) using stellar models which would play the role of 
observed targets. These models included physical phenomena 
not taken into account in the reference models. A total of 13 
targets were constructed, with masses of 0.9 M©, 1.0 M© and 1.1 
M© but to avoid redundancy, we only present 6 of them which 
are representative of the mass and age ranges and of the physical 
effects considered in our study. We tested various effects for 
each mass, namely those that come from microscopic diffusion 
using the approach presented in Thoul et al. (1994) (multiplying 
the atomic diffusion coefficients by a factor given in the last line 
of Table 1^), those caused by a helium abundance mismatch, 
those that result from a metallicity mismatch, as well as those 
which stem from using a different solar heavy elements mixture. 
For the last case, the target was built using the Grevesse-Noels 
1993 (GN93, Grevesse & Noels (1993)) abundances and the 
reference model was built using the Asplund-Grevesse-Sauval- 
Scott 2009 (AGSS 2009, Asplund et al. (2009)) abundances. 

^ These values might seem excessive regarding the reliability of the 
implementation of diffusion. We stress here that our goal was to wit¬ 
ness the impact of significant changes on the results. However, other 
processes or mismatches could alter the ^ gradient and thus be detected 
by the inversion. 


All targets and reference models were built using version 18.15 
of the Code Liegeois d’Evolution Stellaire (CLES) stellar 
evolution code (Scuflaire et al. 2008b) and their oscillations 
frequencies were calculated using the Liege oscillation code 
(Lose, Scuflaire et al. (2008a)). Table 1 summarises the 
properties of the 6 targets presented in this paper. The selection 
of the reference model was based on the fit of the large and 
small separation for 60 modes with n -7-26 and £ - Q — 2 
using a Levenberg-Marquardt minimization code. The use of 
supplementary constraints will be discussed in Sect. 5 whereas 
the effects of the selection of the modes will be discussed in 
Sect. 4. The choice of 60 frequencies is motivated by the number 
of observed frequencies for the system 16Cyg A - 16Cyg B 
by Kepler, which is between 50 and 60 (Vermaetal. 2014). 
The inversions were carried out using the (m, Y) and the (M,ri) 
structural kernels. 

If the inversion of f„ shows that there are differences be¬ 
tween the target and the reference model, then we know that 
the core regions are not properly represented. Whether these 
differences arise from atomic diffusion or helium abundance 
mismatch, the f„ indicator alone could not answer this question'*. 
Therefore, the philosophy we adopt in this paper is simply 
the following: "Is the inversion able to correct mistakes in the 
reference models? If yes, within what accuracy range?". The 
capacity of disentangling between different effects is partially 
illustrated in Sect. 5, but additional indicators are still required 
to provide the best diagnostic possible given a set of frequencies. 

The results are given in Table 2 for the (M,ri) kernels and 
Table 3 for the (u, Y) kernels along with the respective error con¬ 
tributions given according to the developments of Reese et al. 
(2012) and Buldgen etal. (2015). We denote these error con¬ 
tributions: SAvg, Scross, ^Res- These errors contributions being 
defined as follows: 


^-Avg 



(35) 

(36) 


Constraining properly the changes due to multiple additional mixing 
processes with only one structural indicator is of course impossible. 
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Table 1: Characteristics of the 6 targets. The differences between target and reference model are given in boldface. 



Target] 

Targeta 

Target 3 

Target 4 

Targets 

Targets 

Mass (Mo) 

1.0 

1.0 

1.0 

0.9 

1.0 

1.1 

Radius (Rq) 

1.076 

1.159 

1.14 

0.89 

1.193 

1.297 

Age (Gyr) 

8.05 

7.55 

7.06 

6.0 

5.121 

5.135 

Teft (K) 

5597 

5884 

5712 

5329 

6081 

5967 

Zo 

0.015 

0.0122 

0.0122 

0.0122 

0.0122 

0.0122 

To 

0.2457 

0.2485 

0.2485 

0.2485 

0.3078 

0.2485 

Abundances 

AGSS09 

GN93 

AGSS09 

AGSS09 

AGSS09 

AGSS09 

aiviLT 

1.522 

1.522 

1.522 

1.522 

1.522 

1.522 

Diffusion factor 

0 

0 

2 

1.6 

0 

1.6 


if the (m, Ti) couple is used. If one prefers the (u, F), ^cross be¬ 
comes: 


^Cross 



^Cross*^Fr/x. 


(37) 


Finally, erbs is associated to the residual contribution, in 
the sense that it is what remains after one has taken into 
account both Scross and EAvg- The target function and their 
fits are illustrated in Fig. 6 for Target 4 . As we can see from 
Tables 2 and 3, we obtain rather accurate results for all cases. 
This means that the inversion is succesfull and that the regular¬ 
ization process is sufficient for the values (3 - 10“^ and 9 - 10^^. 


We see in the fifth column of Tables 2 and 3 that the aver¬ 
aging kernel fit is usually the dominant error contribution. 
We will see in the next sections how this result changes with 
the modes used or with the quality of the reference model. 
If we analyse the cross-term error contribution, we see that 
it is generally way less important than the averaging kernel 
mismatch error. We also see that despite the high amplitude 
of the cross-term kernels associated with Fi shown in Fig. 6, 
the real error is rather small and often smaller than the error 
associated with the helium cross-term kernel. This is due not 
only to the small variations in Fi between target and reference 
model but also to the oscillatory behaviour of the Fi cross-term 
kernel. In contrast, the cross-term kernel of the (m, Y) kernel has 
a smaller amplitude, but nearly no oscillatory behaviour and 
is more important in the surface regions, where the inversion 
is naturally less robust. Nonetheless, the results of the (m, Y) 
kernels also show some compensation. We also note that they 
tend to have more important residual errors. However, there is 
no clear difference in accuracy between the (m, Fi) and the (u, Y) 
kernels. In the case of Target 2 , we see that although the results 
are slightly improved, the reference value is within the errors 
bars of the inversion results. In an observed case, this would 
mean that the reference model is already very close to the target 
as far as the indicator is concerned. However, we wish to 
point out that it seems rather improbable that the only difference 
between a static model and a real observed star would be in its 
heavy elements mixture. 

Analysing the residual contribution is slightly more diffi¬ 
cult, since it includes any supplementary effect: surface terms, 
non-linear contributions, errors in the equation of state (when 
using kernels related to Y), ... We can see in this study that the 
residual error is rather constricted. This will not be the case 
for example if the parameter 9 is chosen to be very small, or 
if the scaling effect has not been taken into account^. In fact, 

^ In which case one would be searching for a result impossible to ob¬ 
tain. 


the 9 parameter is a regularizing parameter in the sense that it 
does not allow the inversion coefficients to take on extremely 
large values. If that would be the case, the inversion would 
be completely unstable because a slight error on the fit would 
be amplified and would lead to wrong results. This is quickly 
understood knowing that the inversion coefficients are used to 
recombine the frequencies as follows: 


6f^ 



(38) 


with N the total number of observed frequencies. Where this 
equation is subject to the uncertainties in Eq. 17 and Eq. 5 (or, 
respectively Eq. 13), any error will be dramatically amplified by 
the inversion process. Therefore there is no gain in reducing 9 
since at some point, the uncertainties behind the basic equation 
of the inversion process will dominate and lead the method to 
failure. In this case, the inversion problem is not sufficiently reg¬ 
ularized. Such an example is presented and analysed in the fol¬ 
lowing section. 


4. Impact of the type and number of modes on the 
inversion results 

When carrying out inversions on observed data, we are limited 
to the observed modes. Therefore the question on how the 
inversion results depend on the type of modes is of uttermost 
importance. The reason behind this dependence is that different 
frequencies are associated with different eigenfunctions, in other 
words different structural kernels, sensitive to different regions 
of the star. Therefore, the inverse problem will vary for each 
set of modes because the physical information contained in the 
observational data changes. Hence, we studied 4 targets using 
7 sets of modes. As in Sect. 3.3, we wish to avoid redundancy 
and will present our results for one target, namely Target 3 , 
and 5 different frequency sets. As a supplementary test case, 
we defined two target models with the properties of 16CygA 
and 16CygB found in the litterature (Metcalfe et al. 2012; 
White et al. 2013; Vermaetal. 2014). Using these properties, 
we added strong atomic diffusion to the 16CygA model and 
used Z = 0.023 as well as the GN93 mixture for 16CygB. 
In constrast, the reference models used Z = 0.0122 and the 
AGSS09 mixture. The characteristics of these models are also 
summarized in Table 4. We used the set of observed modes 
given by Vermaetal. (2014) and ignored the isolated £ - 3 
modes for which there was no possibility to define a large 
separation. 
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Table 2: Inversion results for the 6 targets using the (m, Ti) kernels. 



^ref 

^ ig^lcm^) (m,Fi) 

^ref 

^ ig^lcm^) 

®Avg 

i(,ri 

^Cross 

wXi 

SRes 

Target; 

4.032 

3.568 + 0.063 

3.532 

4.415 X 10^^ 

-1.684 X 10-4 

9.767 X 10-3 

Target, 

3.434 

3.24 + 0.075 

3.428 

-1.301 X 10-^ 

-4.419 X 10-4 

5.951 X 10-4 

Targetj 

3.562 

3.275 ± 0.067 

3.252 

5.789 X 10-3 

-1.178 X 10-3 

2.277 X 10-3 

Target4 

5.879 

5.621+0.147 

5.536 

1.388 X 10-3 

3.088 X 10-4 

8.062 X 10-4 

Targetj 

2.845 

2.669 + 0.088 

2.630 

-6.498 X 10-4 

-4.493 X 10-3 

2.366 X 10-4 

Targets 

3.205 

3.480 + 0.091 

3.498 

1.496 X 10-3 

-1.11 X 10-3 

7.824 X 10-4 


Table 3; Inversion results for the 6 targets using the (u, Y) kernels. 



^ (g3/cm®) 

^ref 

^ (g^lcm^) (u, Y) 

^ref 

ig^lcm^) 

u,Y 

®Avg 

u,Y 

^Cross 

u,Y 

^Res 

Target; 

4.032 

3.575 + 0.063 

3.532 

8.34 X 10-4 

1.601 X 10-4 

1.1 X 10-3 

Target, 

3.434 

3.423 + 0.075 

3.428 

-1.301 X 10-3 

1.338 X 10-3 

2.127 X 10-3 

Targetj 

3.562 

3.283 + 0.067 

3.252 

5.748 X 10-3 

8.296 X 10-3 

-4.794 X 10-3 

Target4 

5.879 

5.624 + 0.148 

5.536 

1.386 X 10-3 

1.337 X 10-3 

5.448 X 10-4 

Targets 

2.845 

2.675 + 0.089 

2.630 

-5.421 X 10-4 

-8.184 X 10-3 

3.168 X 10-4 

Targets 

3.205 

3.480 + 0.091 

3.498 

1.458 X 10-3 

1.214 X 10-3 

-9.721 X 10-3 



Fig. 6; The left hand-side figure illustrates the averaging kernels for both structural pairs and their respective targets. The red and 
blue curves are nearly identical, so that only the red curve is visible. The right hand-side figure illustrates the cross-term kernels for 

both structural pairs, the target being 0. 


All the sets used for the test cases of this section are sum¬ 
marized in Table 5. The reference model was chosen as in 
Sect. 3.3, using the arithmetic average of the large and small 
frequency separations as constraints for its mass and age. Using 
these test cases, we first analyse whether the inversion results 
depend on the values of the radial order n of the modes, with 
the help of the frequency sets 1, 2 and 3 (see Table 5). Then 
we analyse the importance of the £ - 3 modes for the inversion 
using set 4 and 5. 

The inversion results for all these targets and sets of mode 
are presented in Tables 6 and 7 for the (m, Fi) kernels. A first 
conclusion can be drawn from the results using sets 1, 2 and 3: 
low n are important to ensure an accurate results. Indeed, set 
2 provides much better results than set 1, and even using set 3 
(which is only set 1 extended up to n = 34 for each £) does not 
improve the results any further. This means that modes with 
n> 21 are barely used to fit the target. 

This first result can be understood in different ways. First, 


using mathematical reasoning, we can say that the kernels 
associated with higher n have high amplitudes in the surface 
regions and are therefore not well suited to probing central 
regions. Another way of understanding this problem is that when 
we use modes with high n, we come closer to the asymptotic 
regime and the eigenfunctions are described by the JWKB 
approximation, all having a similar form and therefore not 
providing useful additional information. In that sense, there is 
a clear difference between inverted structural quantities and the 
information deduced from asymptotic relations, which requires 
high n values to be valid, thereby highlighting the usefulness 
of inversions. The question of the importance of the modes 
€ - 3 IS, also quickly answered from the results obtained with 
set 4 and 5. For those test cases, we reach a very good accuracy 
even without n < 9, unlike the previous test case using set 2. 
Moreover, we use even less frequencies than for the 3 first sets. 
This is in fact crucial to determine whether one can apply an 
inversion in an observed case, since a few t - 3 modes can 
change the results and make the inversion successful. 
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Table 4: Characteristics of the 16CygA and 16CygB clones. 



16CygA-clone 

16CygB-clone 

Mass (M 0 ) 

1.11 

1.06 

Radius (Rq) 

1.13 

1.04 

Age (Gyr) 

6.9 

6.7 

Left (K) 

5696 

5772.9 

Zo 

0.024 

0.023 

To 

0.241 

0.242 

Abundances 

AGSS09 

GN93 

aiviLT 

2.0 

2.0 

Diffusion coefficient 

2 

0 


To further illustrate the importance of the octupole modes, 
we use the 16CygA and 16CygB clones and carry out inversions 
for their respective observed frequency sets. First, we use all 
frequencies and reach a reasonable accuracy for both targets. 
In the second test case, we do not use the octupole modes and 
one can observe a drastic change in accuracy. These results are 
illustrated in Table 7, where the notation "Small" (for small 
frequency set) has been added to the lines associated with the 
results obtained without using the ( -?> modes. 

Looking again at the results for 16CygA, we again see 
that although the inversion improves the value of f„, the ref¬ 
erence value lies within the observational error bars of the 
inverted result. The case of the truncated set of frequencies is 
even worse, since the inverted result is less accurate than the 
reference value. We therefore analysed the problem for the full 
frequency set. To do so we carried out different inversions using 
higher values of 0, the results for 6 - 10 are illustrated in 
Table 7. In this case we have smaller error bars, but what is 
reassuring is that the result did not change drastically when we 
changed 0. This means that the problem is properly regularized 
around 0 - 10“^ and 6 - 10 and that we can trust the inversion 
results. Our advice is therefore to always look at the behaviour 
of the solution with the inversion parameters to see if there is 
any sign of compensation or other undesirable behaviour. In fact 
there is no law to select the value of 6 and applying fixed values 
blindly for all asteroseismic observations is probably the best 
way to obtain unreliable results. 

The case of the small frequency set is even more intrigu¬ 
ing since the result improves greatly with 9 - 10^“*. The 
question that arises is whether the problem is not properly 
regularized with 9 - 10“^ or whether we are facing some for¬ 
tuitous compensation effect that leads to very accurate results. 
If we are facing fortuitous compensation, taking 9 slightly 
larger than 10 or increasing p will drastically reduce the 
accuracy since any change in the linear combination will affect 
the compensation. However, if we are facing a regularization 
problem, the accuracy should decrease regularly with the change 
of parameters (since we are slightly reducing the quality of the 
fit with those changes). We emphasize again that one should not 
choose values of the inversion parameters for which any small 
augmentation of the regularization would drastically change 
the result. In this particular case, we were confronted with 
insufficient regularisation and choosing 9 - 10 corrected the 
problem. 


5. Impact of the quality of the forward modelling 
process on the inversion results 

In this section we present various inversion results using dif¬ 
ferent criteria to select the reference model for the inversion. 
The previous results, only using the average large and small fre¬ 
quency separations as constraints for the mass and age of the 
model, are indeed a crude representation of the real capabilities 
of seismic modelling. It is well known that other individual fre¬ 
quency combinations can be used to obtain independent infor¬ 
mation on the core mixing processes and that we should adjust 
more than two parameters to describe the physical processes in 
stellar interiors. To carry out these test cases, we build two tar¬ 
get models including microscopic diffusion. As previously done, 
we did not include this process in the reference models obtained 
by fitting the 56 “observed” frequencies. We used the modes 
Z = 0,n = 12 - 25; f = l,n = 11 - 25; Z = 2,n = 11-26; 
I - ?>,n - 14 - 24. The characteristics of the targets are sum¬ 
marized in Table 8 along with those of the best models ob¬ 
tained through seismic modelling. Table 9 contains information 
on the various constraints and free parameters used for the fit. We 
used various seismic constraints such as the individual large and 
small frequency separations, the individual roi and ro 2 , defined 
as: 

Vn-\,2 - Vnfi 

m =-, 

Vn,l - 

_ Vrifi - 2v„,I + Vn+1,0 
2 (V;i 4-1 0 ^n.o) 

The free parameters were chosen to match what is done when 
trying to fit observations, although the final quality of the fit 
is much higher than what one expects from an observed case, 
as is illutrasted in Fig. 7 for Targetv® and Tab. ?? for Targets, 
where one should note that the arithmetic average of the large 
frequency separations were fitted to within 1 % in addition to the 
individual quantities plotted in Tab. ??. In all these cases the 
inversion improved the value of f„. In some cases, the acoustic 
radius was not well fitted by the forward modelling process 
but the inversion could improve its determination. Table 11 
illustrates the results of combined inversions whereas Table 10 
gives the results for all test cases. One could point out that the 
error contributions are not that different from what was obtained 
using only the average large and small frequency separations as 
constraints. However, our previous analysis of the dependence 
in the degree and radial order of the modes has shown us that 
having low n and I = 3 modes was the best way to ensure 
accuracy. In that sense, the frequency set used for these test 

® We did not present the fit of the individual large frequency separa¬ 
tions for Modeb.i to avoid redundancy with Model 7.2 and Modelv 3 . 


(39) 

(40) 
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Table 5: Sets of modes used to analyse the impact of the number and type of frequencies on the inversion results 



Sell 

Set2 

Set3 

Set4 

Set5 

16CygA 

16CygB 


0 


5-21 

n = 9 - 34 

n = 11 - 24 

n = 11 - 26 

n 

= 12-27 

n = 13 -26 


1 

« = 9 - 28 

5-21 

n = 9 - 34 

« = 11 - 24 

« = 11 - 26 

n 

= 11-27 

n = 13 -26 


2 

n^9-2d, 

5-21 

n = 9 - 34 

n = 11 - 24 

n = 11 - 26 

n 

= 11-24 

n = 12-25 


3 

- 

- 

- 

0 

(N 

1 

o^ 

II 

n = 12-22 

n 

= 15-21 

n = 17 - 24 


Table 6: Inversion results for Target 3 using the (u, Ti) kernels and sets 1 - 5 of Table 5. 



ref 

ref 


wXi 

®Avg 

uXi 

®Cross 

wXi 

Sell 

5.855 

5.700 + 0.161 

5.538 

2.805 X 10-^ 

-4.246 X 10-2 

4.881 X 10-4 

Set2 

5.888 

5.566 + 0.088 

5.538 

4.062 X 10-^ 

1.869 X 10-4 

3.505 X 10-4 

Set3 

5.895 

5.690 + 0.146 

5.538 

2.79 X 10-2 

9.25 X 10-® 

6.8 X 10-4 

Set4 

5.886 

5.570 + 0.110 

5.538 

6.074 X 10-2 

2.859 X 10-4 

-1.893 X 10-2 

Sets 

5.968 

5.630 + 0.105 

5.538 

1.644 X 10-2 

-4.072 X 10-4 

4.714 X 10-4 



Observed frequency fiHz Observed frequency {iHz 


Fig. 7; Results of the fit using the Levenberg-Marquardt algorithm for the first target as a function of the observed frequency; the 
upper panel is associated with Modely i which used the average large frequency separation and the individual roi and dv„i as 
constraints; the central panel is associated with Modelv 2 which used individual Av„_; and §v„j as constraints; the lower panel is 
associated with Model 7.3 which also used individual Av„,; and dv„j as constraints. 


cases is of lower quality and but this has been compensated 
for by the forward modelling process. To further illustrate the 
importance of the model selection, we repeat in Table 10 the 
results obtained by simply ajusting < Av > and < 6v > that we 
denote as Modelv q. We see from the results for Targety that the 
dominant source of error contribution, CAvg is +6 to 10 times 
smaller than what is for the Modely 0 . It is thus clear that using 
the information given by individual frequencies is crucial to 
ensure accurate results in observed cases. 


The necessity of an acoustic radius inversion results from 
two aspects in the selection of the reference models. The first 
one is present in Modely 2 ; the change of a induced during 
the fit had an important impact on the upper regions and 
therefore a change in the acoustic radius was observable. The 
second one is present in Modelg 1 and Model8 2 , where the 
observational constraints were sensitive to core regions, except 
for the arithmetic average of the large frequency separations. 
In this case, the upper regions are thus less constrained and the 
inversion is still necessary. However, we note that in most cases 
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Observed frequency fxHz Observed frequency fj,Hz 


Fig. 8; Results of the fit using the Levenberg-Marquardt algorithm for the second target as a function of the observed frequency 
(colour online); the red dots are associated with Models i which used individual ro 2 as constraints and the purple dots are 
associated with Modelg 2 which used individual ro 2 and roi as constraints. 


the acoustic radius of the reference model was very accurate. 
This is due to the lack of surface effects in the target models. If 
we were to include non-adiabatic computations or differences 
in the convection treatment, differences would be seen in the 
acoustic radii of the targets and reference models but it is 
clear from these test cases that the acoustic radius combined 
with tu alone is not sufficient to disentangle between effects 
of differences in helium abundance and effects of microscopic 
diffusion. We also note that when using the individual ro 2 along 
with the < Av > (the test case of Modelg 1 ), we obtain a very 
good fit of f„ with the reference model. This is a consequence of 
the fact that ro 2 is very sensitive to core regions and therefore 
the core characteristics are well reproduced. However, with the 
acoustic radius inversion, we note that the surface regions are 
not well htted, and if we include roi, as is done in the test case 
of Modelg. 2 , we obtain a better ht of the acoustic radius, but f„ 
is then less accurate. The inversion of the indicator f,, informs 
us that the core regions are not well reproduced in this reference 
model. We can see that Modelg .1 or Model 7.2 and Modelvg 
reproduce better the indicator f„. However, in all these cases 
the acoustic radius was not properly reproduced. Therefore the 
combined inversions indicate that something is wrong with 
the set of free parameters used because we cannot ht properly 
surface and core regions simultaneously, although the ht of the 
seismic constraints with the Levenberg-Marquardt algorithm is 
excellent in all these cases. 

We also mention that in all test cases carried out here, we 
did not consider the hrst age indicator t from Buldgen et al. 
(2015). In fact, the indicator t could also provide accurate results 
for Targetv. However, its inaccuracy for older models has been 
observed during this entire study and we recommend to limit its 
use to young stars^, where it can provide valuable information 
provided the kernels are well optimised. 


6. Conclusion 

In this article, we have presented a new approach to constrain 
mixing processes in stellar cores using the SOLA inversion 

’’ We chose not to present these results to focus our study on the 
indicator. 


technique. We used the framework presented in Buldgen et al. 
(2015) to develop an integrated quantity, denoted f„, sensitive to 
the effects of stellar evolution and to the impact of additional 
mixing processes or mismatches in the chemical composition of 
the core. We based our choice solely on structural effects and 
considerations about the variational principle and the ability of 
the kernels to ht their targets. 

The derivation of this new quantity was made possible by 
the use of the approach of Masters (1979) to derive new struc¬ 
tural kernels in the context of asteroseismology by solving an 
ordinary differential equation. We discussed the problem of the 
intrinsic scaling effect presented in Basu (2003) and discussed 
how it could affect the indicator f„. We tested its sensitivity to 
various physical changes between the target and the reference 
model and demonstrated that SOLA inversions are able to 
signihcantly improve the accuracy with which f„ is determined, 
thereby indicating whether there is a problem in the core regions 
of the reference model. 

We also analysed the importance of the number and type 
of modes in the observational data and concluded that the 
accuracy of an inversion of the indicator f„ increased with 
multiple values of the degree, {, and low values of the radial 
order, n. In that sense, we emphasize that the observation of 
{ - 3 modes is important for the inversion of the indicator 
since it can improve the accuracy without the need of low n 
modes. Such modes are difficult to observe. Indeed, only a few 
octupole modes have been detected for around 15% of solar-like 
stars with Kepler. The use of other observational facilities, such 
as the SONG network (Grundahl et al. 2007), might help us 
obtain richer oscillation spectra as far as the octupole modes are 
concerned. The test cases for the 16CygA and 16CygB clones 
demonstrated that our method was applicable to the current 
observational data and one could still carry out an inversion of 
without these modes. However, it is clear that this method will 
only be applicable to the best observational cases with Kepler, 
Plato or SONG. 

We also analysed the impact of the selection of the refer¬ 
ence model on the inversion results and concluded that using 
individual frequency combinations is far more efficient in terms 
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of accuracy and stability of the inversion results. However, we 
also noticed in supplementary tests that there was what could 
be called a “resolution limit” for the inversion which depends 
on the magnitude of the differences in the physics between the 
reference model and its targets but also on the weight given 
to the core in the selection of the reference model (see for 
example the case of Models i and the associated discussion). 
This leads to the conclusion that supplementary independent 
integrated quantities should be derived to help us disentangle 
between various physical effects and improve our sensitivity to 
the physics of stellar interiors. Nevertheless, the test cases of 
Sect. 5 showed that the SOLA method was much more sensitive 
than the forward modelling process used to select the reference 
model (here a Levenberg-Marquardt algorithm) and could 
indicate whether the set of free parameters used to describe the 
model is adequate. 
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Table 7: Inversion results for the 16CygA and 16CygB clones using the (u, Ti) kernels. 






uXi 

®Avg 

uXi 

®Cross 

wXi 

16CygA (Full, 9 = lO^'’) 

2.965 

2.891 +0.083 

2.885 

2.641 X 10-^ 

-1.780 X 10-3 

1.442 X 10-3 

16CygA (Full, 0 = 10-"^) 

2.965 

2.872 ± 0.036 

2.885 

-4.033 X 10-3 

-9.487 X 10-'' 

7.149 X 10-4 

16CygA (Small, 0 = lO^^) 

2.965 

2.971 +0.083 

2.885 

3.117 X 10-3 

-2.577 X 10-3 

5.000 X 10-4 

16CygA (Small, 0 = lO^^) 

2.965 

2.906 + 0.031 

2.885 

8.240 X 10-3 

-1.778 X 10-3 

7.000 X 10-4 

16CygB (Full) 

4.540 

4.007 + 0.095 

3.783 

4.547 X 10-3 

-2.407 X lO-** 

2.277 X 10-3 

16CygB (Small) 

4.540 

4.295 + 0.113 

3.783 

1.093 X 10-' 

-1.715 X 10-3 

1.138 X 10-3 


Table 8 ; Characteristics of Targetv, Targets and of the models obtained with the Levenberg-Marquardt algorithm. 



Targetv 

Model? 1 

Model?.? 

Model?.? 

1 Targets 



Mass (M 0 ) 

0.9 

0.933 




1.009 

1.029 

Radius (Rq) 

0.908 

0.919 


0.926 

1.17 

1.18 

1.19 

Age (Gyr) 


3.34 

3.69 

3.45 

4.168 

4.322 

4.489 

Teff (K) 

5659 

5701 

5488 

5713 


5985 

5966 

Z 


0.0122 


0.0105 


0.0185 

0.0181 

To 


0.274 


0.243 


0.323 

0.305 

Abundances 


AGSS09 


AGSS09 

BnrSSnS 

AGSS09 

AGSS09 

aM+T 

1.522 

1.522 

1.297 

1.522 

1.522 

1.522 

1.522 

Diffusion coefficient 

1.6 

0 


0 

1.6 

0 

0 


Table 9: Constraints and free parameters used for the Levenberg-Marquardt fit. 



Constraints 

Parameters 

Model?o 

< Av„./(v) > + < 6 Vn,l(v) > 

Mass + Age 

Model? 1 

Av„,;(v) + §Vn,i(v) + roi(v) 

Mass + Age + To 

Model? 2 

Av„./(v) + 6 Vn,l{v) 

Mass + Age + yo+ ^mlt 

Model? 3 

Av„.;(v) + 6 Vn,l(v) 

Mass + Age + Yq + Zq 

Models. 1 

< Av> + ro2(v) 

Mass + Age + Yq + Zq 

Models.2 

< Av > + ro2(v) + roi(v) 

Mass + Age + Yq + Zq 


Table 10: Inversion results for various fits with the Levenberg-Marquardt algorithm. 






uXi 

®Avg 

uX\ 

^Cross 

wXi 

^Res 


7.626 

7.000 + 0.24 

6.703 

4.952 X 10-4 

2.94 X 10-4 

-3.389 X 10-4 

Model? 1 

6.562 

6.657 + 0.122 

6.703 

-1.4 X 10-^ 

-2.407 X 10-4 

-2.777 X 10-4 

Model? 2 

6.393 

6.669 + 0.123 

6.703 

-6.597 X 10-3 

-1.292 X 10-4 

-2.478 X 10-3 

Model ?.3 

6.467 

6.667 + 0.121 

6.703 

-2.571 X 10-3 

-3.322 X 10-4 

-5.661 X 10-4 

Models. 1 

3.450 

3.651 +0.092 

3.568 

2.458 X 10-3 

-2.302 X 10-4 

-1.745 X 10-3 

Models .2 

3337 

3.637 + 0.096 

3.568 

2.089 X 10-3 

-6.553 X 10-4 

-1.568 X 10-3 


Table 11: Combined (t, f„) inversion results for various fits with the Levenberg-Marquardt algorithm. 






Tref(s) 


Tobs(s) 

Model? 2 

6.393 

6.657 + 0.092 

6.703 

3230 

3223 + 0.028 

3222 

Models. 1 

3.450 

3.651 +0.096 

3.568 

4509 

4450 + 0.028 

4442 

Models .2 

3.337 

3.637 + 0.120 

3.568 

4517 

4448 + 0.017 

4442 
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